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Abstract 

We propose a more direct approach to constructing differential 
operators that preserve polynomial subspaces than the one based on 
considering elements of the enveloping algebra of st(2). This approach 
is used here to construct new exactly solvable and quasi-exactly solv- 
able quantum Hamiltonians on the line which are not Lie-algebraic. It 
is also applied to generate potentials with multiple algebraic sectors. 
We discuss two illustrative examples of these two applications: an 
interesting generalization of the Lame potential which posses four al- 
gebraic sectors, and a quasi-exactly solvable deformation of the Morse 
potential which is not Lie-algebraic. 

1 Introduction 

Our purpose in this paper is to show that the property of quasi-exact solvabil- 
ity for Schrodinger operators admits an effective formulation which goes be- 
yond the Lie algebraic context of hidden symmetry algebras, and to show ex- 
plicitly that several classes of physically relevant potentials are quasi-exactly 
solvable in this more general sense. These include novel generalizations of 
the Lame and Morse potentials. 
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To put our results in context, we begin by briefly recalling the general 
definition of quasi-exact solvability, as well as the Lie algebraic formulation 
which applies to some particular classes of invariant subspaces. 

Consider a one-dimensional Schrodinger operator 

H = -d xx + U(x), (1) 

which is assumed to be essentially self-adjoint on a Hilbert space H of wave 
functions. We say that H is quasi-exactly solvable (QES) [16, 19] if it leaves 
invariant a non-trivial finite-dimensional subspace M. C 7i 

HM C M 

where M is 

M = span{0!, . . . , n }, <f>i E H. 

In this case the restriction of H to Ai is represented by a finite dimensional 
matrix, and a finite portion of the spectrum of H can be obtained alge- 
braically by diagonalizing this matrix. If H actually preserves a complete 
flag of subspaces 

Mi c M 2 c • • • c M n c ... 

then it is said to be exactly solvable (ES). Given a certain potential U(x) it is 
in general very difficult to determine whether H preserves some non-trivial 
invariant subspace At. One generally proceeds in the opposite direction, 
and starting from a known finite-dimensional subspace At one attempts to 
determine all Hamiltonians H which are known a priorito leave M. invariant. 
In most applications M. is assumed to be of the form 

M n = n(x)V n (z(x)), (2) 
V n ={l,z,z 2 ,...,z n ), (3) 

where \i is some smooth non-vanishing multiplier, so that the first problem is 
to construct second-order differential operators which leave the polynomial 
module V n invariant. In this context, Burnside's theorem ensures that every 
operator which leaves V n invariant is an element of the enveloping algebra of 
the Lie algebra sl(2) in the following standard representation 

Jt = z 2 d z -nz, J° n = zd z --, .r n =d z . (4) 
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This approach is known as the Lie algebraic approach and sl(2) is said to be 
a hidden symmetry algebra of the Hamiltonian. The Lie-algebraic approach 
has been extended to higher dimensional problems: some examples of Lie- 
algebraic Hamiltonians in two dimensions are given in [13], while the quantum 
iV-body Calogero-Moser model has been shown to have hidden symmetry 
algebra sl(N + 1) [20], and many-body problems with elliptic interaction 
were shown to be Lie-algebraic in [8]. 

In the last twenty five years, most of the papers dealing with quasi-exact 
solvability have made use of the Lie-algebraic approach, using irreducible 
representations of finite-dimensional Lie algebras by first-order differential 
operators as building blocks for constructing differential operators with in- 
variant subspaces. Despite the success of this method, there remain some 
natural questions which are not amenable to this Lie-algebraic approach. 
For example, one would like to be able to construct the most general differ- 
ential operator that will map a polynomial module V n into a proper subspace 
V n -k C V n or to know what are the operators that preserve monomial sub- 
spaces, that is subspaces of V n which admit a monomial basis, but where all 
degrees up to n need not be present (see definition (JHJ) in Section^. It is the 
latter question that we consider in Section 2 of this paper, where we obtain 
a general structure theorem for the space of linear differential operators pre- 
serving a general monomial subspace. Much of the contents of this Section 
can already be found in a paper by Post and Turbiner in [18]. Sections 3 and 
4 form the heart of our paper. Section 3 consists in an explicit application 
of the results of Section 2 to Hamiltonians with multiple algebraic sectors 
viewed from the perspective of monomial subspaces/ This Section also in- 
cludes an intriguing generalization of the Lame equation which is shown to 
admit four algebraic sectors instead of the double algebraic sector admitted 
by the classical Lame equation. In Section 4 we discuss a modification of 
the Morse potential which is shown to preserve a proper monomial subspace 
of V n . This class of potentials cannot be obtained in the Lie-algebraic ap- 
proach and therefore they are not present in the classification of Lie-algebraic 
potentials on the line [11]. 
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2 The direct approach to invariant monomial 
subspaces 

In this Section we shall introduce the concept of a monomial subspace and 
we will characterize all differential operators of order r that preserve a given 
monomial subspace. This Section follows essentially the algebraic approach 
introduced in [18]. 

For a finite subset / C Z let us define a monomial subspace Vi as 

Vi = span{z fc : k e I}, (5) 

and let N = \I\ be the dimension oiVi- A non-trivial differential operator 
Tf of order r will be said to be of degree d G Z if and only if 

T r d [x j ] = Cj x j+d Vj'eZ. (6) 

Let us denote by D r (Vi) the space of differential operators of degree r with 
analytic coefficients that preserve Vi- 

Proposition 1 Every linear transformation of Vi can be represented by an 
element of D^^iiVi). 

Proof. A basis of the dual space is given by the following (N — l)-th order 
operators 

^=^nfe^. (7) 

jei J 

which clearly satisfy L\\z^\ = 5 J k , Vj, k G /. Therefore, every element of 
End('Pj) can be represented by an operator of the form 

j,kei 

QED 

In fact, these operator duals can be constructed for more general sub- 
spaces of functions (not just polynomials) [15], but the order of the differen- 
tial operators grows with the dimension of the space. We are interested in 
constructing differential operators that preserve Vi of order r < N — 1. Let 
us denote by Df(Vi) C D r (Vi) the space of operators of order r and degree 
d that preserve Vi. In order to characterize Q r (Vi) we can concentrate in 
operators having a fixed degree, since the following Proposition holds. 



4 



Proposition 2 

(8) 

d 

proof Any operator T 6 ® r (Vi) can be decomposed uniquely as a sum 
T = ^ d T^ d ' of operators of fixed degree. Then it suffices to prove that if T 
preserves Vi then each T"W also preserves Pj. Since TP/ C Vj, then 

= ^2a jk z J , kel, 




a d +k,kZ d+k , if d + k el; 

otherwise. 

which proves that T^Vi C P/.QED 
Let us define the primitive index sets 

I(a,b,c) = {a + kc: < k < b} (9) 

and the primitive monomial subspaces V a ,b,c as 

V aM = V I{ a, b ,c) = span{z a+fcc :0<k<b} (fO) 

where a is the smallest exponent, c is the size of the gap and b + 1 is the 
number of elements. In fact, 

V aAc = z a P b (z c ), 

where V n {z) = (1, z, . . . , z n ) denotes the full polynomial module (no gaps). 
Every monomial subspace can be trivially decomposed as a direct sum of 
singleton primitive monomial subspaces 



0, 
a£l 

which corresponds to representing / as the union of iV singletons. The dif- 
ferential operators (J7J) of order N — 1 are constructed to annihilate N — 1 
different exponents. However, it is precisely when the index set / can be 
expressed as the union of a smaller number of primitive index sets ©, that 
lower order operators preserving Vi can exist. This consideration leads to 
the following key theorem: 
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Theorem 1 If the set I can be expressed as the union of a finite number of 
primitive index sets of the form 

r 

I = \Jl(ai,bi,d), (11) 

i=i 

then there exist operators Tr ±d ^ of order r and degree ±d that preserve Vi- 
Conversely, if is an operator of order r < N and degree d > that 
preserves a given monomial subspace Vi, then there exists a Tj d ' which also 
preserves V i and a vector-space decomposition ofVi into primitive monomial 
subspaces 

r 
i=l 

Proof. Suppose that (jllj) holds and consider the following r-th order 
operators 

r 

T y d ) = z- d l[(zd z -a l ) (12) 

i=l 
r 

= z d H(zd z -(a i + db i )) (13) 

i=l 

The lowering operator maps z ai+kd to z ai+ ( k ~ 1 > d and annihilates z ai for all i = 
1, . . . ,r. Similarly the raising operator maps z ai+kd to z ai+( - k+1 ^ d annihilates 
z ai+dbi £ or ^ jt j s c i ear then that 

T r (±d) P/ c Pi. 

To prove the converse we assume that is an operator of degree d > 
and order r that preserves a given Vi. The operator z~ d T^ has degree 0, 
and hence has a unique factorization of the form 

r 

z- d T^ = cl[(zd z -k i ), 

i=i 

up to a multiplicative constant c. For each i, let bi to be the largest possible 
natural number such that I(ai, bi, d) C /, where = hi — db{. 
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If (jllj) did not hold, suppose that a is the largest element of I that is not 
in any of the I(di,bi,d). Since a 7^ k{ for all i, T^ d ' does not annihilate z a 
and by assumption a + d G /. But since a was the largest element in / which 
does not belong to any of the I(cii,bi,d), we must have a + d G I(di,bi,d) 
for some i, but then by the way 6, was chosen, a G / as well, which is a 
contradiction. QED 

2.1 Second order operators 

The case of second order operators is particularly important for the appli- 
cations in quantum mechanics, and we shall concentrate to their study in 
this Section. From Theorem ^ we know that if a given monomial subspace 
Vi is invariant under T 2 d \ it is also invariant under T 2 , and of course it 
always happens that T 2 ^ G D 2 (Vi) for any monomial subspace Vj. The 
natural question is now to explore whether a given monomial subspace can 
be preserved by several operators of different degrees. 

Theorem 2 Suppose that a monomial subspace Vi of dimension \I\ > 3 is 
invariant under two second order operators and of degrees d\,di G 
Z + with d\ > d 2 coprimes. Then d\ = 2, d 2 = 1 and Vi must be one of the 
following 3 modules: 

l) V n = (l,Z,Z 2 ,...,Z n ) 
%%) V n = (l,Z, Z \...,Z n -\z n ) 

III) V n = {l } z\... } z n ) 

Proof. Let / = {ii, . . . ,i^} with i 1 < i 2 < ■ ■ ■ < in-i < i>N- Both T2 
and T^ 6,2 ' 1 must annihilate z %N , so they must be of the form 

T 2 {dl) = z d ^(zd z - i N )(zd z - ai ) 
T 2 {d2) = z d *(zd z -i N )(zd z -a 2 ) 

with ai and a 2 still to be specified. The difference ix — iN-i must be equal to 
either d\ or d 2 because otherwise both T 2 and T 2 d2 ^ must annihilate z lN ~ 1 
too, fixing ai = a 2 = i N _ x . But then either T 2 {dl) [z iN ~ 2 ] or T 2 {d2 \z iN - 2 } will be 
outside Vi- So two possibilities remain 
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1. i N - i N _i = di 

In this case must annihilate z tN ~ 1 so a 2 = ijv-i, but this forces 
tN-i — in-2 = d 2 since otherwise T^f^- 2 ] ^ Vj. It follows that 
T 2 {dl) [z iN - 2 ] = fixing ai = ijv_ 2 . But now neither T 2 (d2) nor T 2 (dl) can 
annihilate z JJV ~ 3 , so in-2 — = ^2 and ijv-i — 2jv-3 = di- This 
implies di = 2d 2 and since they are coprimes d± = 2 and d 2 = 1. The 
preserved module is then 7\. 

2. i^v — ^n-i = d 2 

In this case T^ 1 ^^ -1 ] = so a\ = In-i- The difference In-i — i>N-2 
cannot be d\ because otherwise would need to annihilate z lN ~ 2 , 
and then either T$ dl) [z itf - 3 ] or T 2 {d2) [z iN - 3 ] £ Vj. So i^-i - i N _ 2 = d 2 
but then T 2 ^' l \z lN -' 2 \ oc z lN and this can only happen if d\ = 2d 2 . 
Iterating these arguments it is easy to see that all differences between 
successive indices must be equal to d 2 except for i 2 — i\ which could be 
either d\ (which leads to V n ) or d 2 , which leads the ordinary polynomial 
module P„.QED 

Corollary 1 No monomial subspace Vi is invariant under three second order 
operators of relatively prime degrees di, d 2 , d 3 e Z + . 

Proof. It is clear that a third operator T 2 d3 ^ with an odd d 3 > 2 will not 
preserve any of the modules of the previous Theorem. 

Corollary 2 The spaces of second order differential operators that leave V n , 
V n and V n invariant have the following dimensions: 

dim!D 2 (7\) = 9, dimD 2 (P„) = dimD 2 (P n ) = 7. 

Proof. For the case of V n , the dimensions of the corresponding operator 
spaces D 2 (V n ) are: 

dimS)°(P n ) = 3, dimDj 1 ^) =2, dim2)± 2 (P n ) = 1. (14) 



8 



This follows from noting that the operators in question can be expressed in 
the form 



T (o) 
T (i) 
T (-i) 

T (-2) 
T (2) 



p{zd z ), p{x) e V 2 [x) 
z{a\zd z — bi)(zd z — n) 
z~ l (a 2 zd z - b 2 )(zd z ), 



v-2, 



zd z - l)(zd z 



z 2 (zd z - n)(zd z - (n- 1)). 



The situation for the case of V n is similar, the difference being that the 
dimension of 'S^vn) is only one since the corresponding operators must 
have the form 

T (1) = a lZ (zd z -0)(zd z -n), 
T ( - x) = a 2 z'\zd z -0)(zd z -2). 

Thus the dimension of D 2 {V n ) is seven. QEDIn order to apply this con- 
struction to Quantum Mechanics, we need to relate an operator T e 2) 2 (P/) 
with a Schrodinger operator (JIJ. This is always possible in one dimension by 
virtue of the following Lemma: 

Lemma 1 Every second order operator T of the form 

T = P(z)d zz + (Q(z) + l -P'{z))d z + R(z), (15) 

with P(z) > can be transformed into Schrodinger form (JIJ by a gauge 
transformation and a change of variables: 



H(x) = -fi(z)T(z)^ L (z 



where 



x 



fi{z) 



exp 



Q(C) 



2P(C) 



The resulting potential U(x) in the physical variable x is given by 



(16) 

(17) 
(18) 

(19) 
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3 Quantum Hamiltonians with multiple alge- 
braic sectors 



We are familiar with the fact that Schrodinger operators with even poten- 
tials will have eigenf unctions with well defined parity, and therefore the odd 
sector can be treated independently from the even sector. Let us look at the 
phenomenon of multiple algebraic sectors from the point of view of monomial 
subspaces. There are two mechanisms that can produce a multiple algebraic 
sector. The first one comes from the fact that a second order operator can 
preserve two different polynomial modules. From Theorem \I\we know that 

can be decomposed as 

V I = P ai , m ,d + 'Pai,n2,d- (2°) 

If in addition D^2 2 \Vi) = for every other d 2 ^ d,d 2 G Z + then 

3>£*(P/)c!D£Vai,»i,d), (21) 
sf(7>i)c£f(7V n2)d ), (22) 

and both spaces are invariant separately. We can take a± = without loss 
of generality and let p = a 2 /d,p G R. In the new variable y = z d the two 
invariant polynomial modules are 

Mx=P ni (y), M 2 = y p V n M, (23) 
and each of them is independently preserved. 

Proposition 3 The most general second order differential operator that pre- 
serves both Mi and M. 2 in (|23*j) is a linear combination of the following five 
operators 

T 2 (+1) =y(yd y -n 1 )(yd y -(p + n 2 )) (24) 

= y 3 d yy + (1 - p - ni - n 2 )y 2 d y + n x (n 2 + p)y, 
if =y*d yy , (25) 

Tt l) =y-\ydy-0)(yd y -p) (26) 

= ydyy+(l-p)dy, 

Tf =yd y , (27) 
T (0) = 1. (28) 
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These operators form a sub-family of the well known nine parameter family 
of Lie-algebraic second order operators that would preserve the polynomial 
module M.\. 

A second source of multiple algebraic sectors comes from noting that 
in Lemma ^ the correspondence between the algebraic operator T and the 
Schrodinger operator H is at least two to one: 

Proposition 4 The operators T\ and T 2 given by 

Ti =P{z)d zz +(Q{z)+ l -P'{z))d z + R{z) 1 (29) 
T 2 = P(z)d zz + ( - Q(z) + \p'{z))d z + R(z) - Q'(z) + l -P"(z) (30) 

are both equivalent under the transformation - (|18|) to the same Schrodinger 
operator H = —d xx + U(x) with potential U(x) given by (fT9*j) . 

Therefore if both 7\ and T 2 preserve different finite dimensional subspaces, 
then the Schrodinger operator will have a double algebraic sector. 

The two mechanisms described above can be combined to yield potentials 
with quadruple algebraic sectors. This will be the case if operators T\ and T 2 
in Prop osit ion 0] are both a linear combination of the operators of Proposition 
El These ideas lead to a generalization of the Lame potential that posses a 
quadruple algebraic sector. 



3.1 The generalized Lame potential 

Let us consider the following linear combination of the operators (|24j) - (J28j) : 

T = 4fcT 2 (+1) + 4(k + 1)T 2 (0) + 4T 2 M) + 46 1 T 1 (0) (31) 

which can be written in the form (fl3|) with 

P(z) =4z(z + l)(kz + l), (32) 
Q(z) = -2k{l + 2m + 2n 2 + 2p)z 2 + 4(6 X - A; - l)z + 2 - 4p, (33) 
R(z) =Akn 1 (n 2 + p)z. (34) 

The cubic P(z) determines the change of variables (|17|) and we must consider 
only the interval in z such that P{z) > for the change of variables to be 
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admissible. If we take k < and let z G (0, —1/k) the change of variables 
determined by (fT7j) is 

z = k~ l sn 2 {^/^kx\k- 1 ), (35) 

which maps the real line to the interval (0,-1/k). From equation (|T9~)) in 
Lema^we see that the condition for the potential to be nonsingular is that 
Q(z)(P'(z) — Q(z)) should divide P(z), and since z G (0, —1/k) this implies 
that 

Q(0) = or Q(0)=P'(0), 

and 

Q(-l/k) = or Q{-l/k) = P'{-l/k). 

The first of these conditions is satisfied when p = ±1/2 and since we can as- 
sume, without loss of generality, that is the smallest power in the monomial 
subspace, then we must take p = 1/2. The second set of conditions amounts 
to choosing b\ as 

bi = k — rii — n 2 , (36) 
h = 2k-l-n x -n 2 . (37) 

Therefore we can write the potential (fTTIJ) after performing the change of 
variables as 

V(z) = ml(l + l)sn 2 (a;|m) - j{j + l)sn 2 (Vm- lxl^) (38) 

where m = 1 — k and j, I e Z are related to n\ and n 2 as 

I — m +n 2 + 1, j = ni — n 2 — 1 (39) 
/ = m + n 2 + 2, j = n 2 - n x (40) 

depending on whether we use fJBljj) or ([3T)l. Note that the particular cases 
j = and j = — 1 correspond to the Lame equation, and in those cases the 
two polynomial modules in (}2*3*j) match with no half-integer gaps. However, 
for greater values of j the potential departs from Lame and develops an 
interesting shape, as can be seen in Figure ^ Depending on whether we use 
(HSU) or (|3*7j) in the formula for the gauge factor (fTSj) we obtain 

A t 1 = (1 + z )- l/2 = cn l (x\m), (41) 
^2 = (1 + z)- l/2 (l + kz) 1/2 = cn'(x|m)cn( v /Z A?x| T ^), (42) 
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Figure 1: Generalized Lame potential (jHEJ) for m = 11/10, / = 9 and j = 6 



To summarize the results in this section, we have proved that when j and / 
are integers of the same parity, the generalized Lame potential ()38|) has the 
following quadruple algebraic sector: 



where 



Aii V ni (z) 
1+3 



I 



2 

1-3 



-1, 



n 2 



n 2 



-- HxV^iz), 
3 1 



2 

l + 3 



(43) 
(44) 



(45) 



and z is given by t)35j) . The total number of algebraic eigenfunctions is the 
sum of the dimensions of the four algebraic sectors, which equals 2l + 1. It is 
quite remarkable that this number does not depend on j. In fact, when j = 
and I G N, these eigenfunctions are the well known 21 + 1 Lame polynomials. 
The period of the potential is 2K where 



K 



<:W 



yl — m? sin 2 9 



(46) 



is the elliptic integral of the first kind. There are / + 1 eigenfunctions which 
belong to the sectors Aixl have the same period as the potential, while the I 

eigenfuntions which belong to the sectors M.\ 2 have period 4K. A detailed 
study of the properties of these polynomial eigenfunctions shall be deferred 
to a forthcoming publication. 



13 



Let us note as a conclusion that the generalized Lame potential could be 
applied as realistic model of a one- dimensional crystal, which generalizes the 
models introduced by Alhassid et al. in [1]. The generalized Lame poten- 
tial would describe a one dimensional crystal with constant lattice spacing 
a = 2K and two atoms of different atomic number. By changing the values 
of j and / the relative strength of the potential wells can be modified. The 
eigenvalues of the potential can be calculated algebraically and they describe 
the edges of the allowed energy bands of the crystal. Other possible applica- 
tions of this potential include the quantum fluctuations of the inflaton field 
in certain cosmological models [6,14]. 

4 New Quasi-exactly solvable hamiltonians which 
are not Lie-algebraic 

In this section we will apply a similar construction to obtain quasi-exactly 
solvable potentials that cannot be obtained as a quadratic polynomial in the 
generators of si(2). In order to construct second order operators that are not 
Lie-algebraic we should focus our attention in cases ii) and Hi) of Theorem 
121 the so called exceptional modules, since the full polynomial module V n is 
the carrier space for an irreducible representation of st(2), and we have from 
Burnside's theorem that any operator which leaves V n invariant must be a 
polynomial in the generators of st(2). Let us note that the two exceptional 
modules are projectively equivalent, since 




(47) 



and we can therefore restrict our attention to one of them. 



Proposition 5 A differential operator of order two or less preserves 




(48) 
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if and only if it is a linear combination of the following seven operators: 



T (+2) 
1 2 


= z 4 d zz 


+ 2(l-n)z 3 d z + n(n-l)z 2 


T (+l) 
1 2 


= z 3 d zz 


+ (l-n)z 2 d z , 


T (0) 
1 2 


= z 2 d zz 


j 


1 2 


= zd zz - 


- d z , 


T (-2) 
1 2 


= d zz - 


2z~ l d z , 


T (o) 


= zd z . 




T (o) 


= 1. 





Proof. The sufficiency follows from the fact that each of the operators above 
preserves P n , which can be verified by direct computation. The necessity is 
a direct consequence of Corollary |2J and the fact that the above operators 
are linearly independent. 

Consider the action of SL(2, R) on the (projective) line according to the 
linear fractional, or Mobius transformations 

w ^z = aw + (3 A = a<J-/?7 = l. (49) 
jw + o 

There is an induced action of SL{2) on the space V n of polynomials of degree 
at most n, given by 

P(z) i — > P(w) = («fw + 5) n P ( aw + P \ . (50) 



•yw + 5 

This irreducible multiplier representation of SL(2) is isomorphic to the stan- 
dard representation on homogeneous polynomials of degree n in two variables. 
Since V n is invariant under this transformation, it was used in [13] to classify 
the second order operators that preserve Vn into equivalence classes, each 
of them labeled by a simple canonical form. On the contrary, monomial 
subspaces (and in particular the exceptional module V n ) are not invariant 
under the transformation (|5U|). This is clear since V n C V n is a hyperplane 
of co-dimension 1 in V n which moves under the projective SL(2) action. 
In particular, under the transformation (joTIj) the exceptional module V n (z) 
transforms into 

V n (z) !-»• M n (w) C V n (w), 

M n {w) = span { (jw + 5) n ~ k {aw + (3) k \ k = 0, 2, . . . , n) . 
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It will be convenient to introduce a basis independent characterization of the 
space A4 n . 

Proposition 6 //aj^O a polynomial P(w) of degree n is in M. n C V n if 
and only if 

P'(-P/a) - nar/P(-P/a) = 0. (52) 

Proof. Since the condition is linear, it suffices to prove that it is satisfied 
by all the elements in the basis (|51|) of A4 n (w). It is clear that all elements 
with k > 2 satisfy P'(—f3/a) = P(—(3/a) = while the element with k = 
satisfies precisely condition ()52|). In fact, condition ()52|) is the transformation 
of -P'(O) = (which defines V n ) under the transformation ()50p.QED 

We would now want to know which operators preserve the transformed 
monomial subspace M. n . Of course if T(z)V n (z) C V n (z) then 

f(w) = {jw + 5) n -T (?^±A\ . ( JW + 5 )-n 

V 7U> + o J 

will satisfy T{w)M. n {w) C A4 n (w). Transforming the operators of Proposi- 
tion Sunder this prescription will lead to a basis of D 2 (M-n)- However, a 
simpler basis can be written in the form expressed by the following proposi- 
tion: 

Proposition 7 A basis of f £>2{M. n ) is given by the following seven operators: 



Jli 


» 


= (aw 


+ P) 4 d ww + 2a(l - n)(aw + f3fd w + n(n - l)a 2 (aw + f3) 2 


J 2 ( 


» 


= (aw 


+ Pfd ww + a(l - n)(aw + f3fd w , 


J 3 I 


» 


= (aw 


-1- /3) WW , 




» 


= (aw 


+ P)d ww - a(l + nj(aw + 0))d w , 


J 5 ( 




= d 


2a 2na 2 7 
a .d w 2na-fd w + 
(aw + p) (aw + p) 


J 6 I 


» 


= (aw 


+ P) 2 d w - na(aw + f3) + ^— ' + ^ d w , 

7 


J7I 


» 


= 1. 



(53) 

Proof. According to Proposition |Hl a different basis for M. n (w) is given by 
{y + —,y 2 ,...,y n \y = aw + p}, (54) 
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and a straightforward computation will verify that the set of operators 

y% y + 2(1 - n)y 3 d y + n{n - l)y 2 , 
y 3 d yy + (1 - n)y 2 d y , 



y 9yy , 

ydy y - (1 + rvyy)d y , 

dyy - 2y^dy ~ l^dy + , 

y 2 d y -ny+ -d y , 

7 

preserve the vector space spanned by (J33J). Making the appropriate substi- 
tutions yields directly (f53|).QED 

4.1 An example: the modified QES Morse potential 

Consider the following linear combination of the operators f!53|) : 

T(z) = 8J 3 (z) - 8V2J 4 {z) + U 5 (z) - 16J 6 (z). (55) 

If we choose a = (3 = 5 = l/\/2 and 7 = —l/y/2 the above operator can be 
written in the form ()15|) with 

P{z) = 4z 2 } (56) 
Q(z) =-4^2z 2 + (n + l)z-4 + j^j (57) 

An 

R{z) = 8n{l + z)- — — . (58) 

\ ~\~ z 

The change of variables determined by (|17j) is z = e 2x and the potential is 
defined and free of singularities for all values of x. In particular, the potential 
form determined by (fT^j) with (|5T j) - (jo^j) is 

Un(x) = i ^cosh4x - ncosh2x + (1 ^ &2x)2 - ( 59 ) 

where an additive constant has been dropped. If the last two terms were 
absent, this potential would be a QES deformation of the Morse potential, 
which belongs to the classification of Lie-algebraic potentials performed in 
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Figure 2: The deformed Morse potential U n (x) in for n = 2, 10 and 15. 



[11]. A plot of this potential form for different values of n can be found in 
Figure El When n is a non-zero integer, n eigenfunctions of the potential 
()59|) can be calculated algebraically. They have the following form: 

ipk-i{x) = n{x)p h (e 2x ), k = l,...,n (60) 

where fj,(x) is determined by (|TH|) with (J51)j) - (J55j) to be 

-(e- 2x +e 2x +(n-l)x) 

M*) = TT^ (61) 

and Pk(z) is one of the n polynomial eigenfunctions that the operator (|B3|) 
has in the space 

n — 2 

M n {z) = span{^ + — — , (z + l) 2 , ...,{z + l) n }. (62) 

To be more specific, let us calculate a few eigenfunctions explicitly. If we let 
n = 3, the action of T(z) with respect to the basis (|6^)l of M.zi. z ) is given 
by: 

/ -12 6V2 \ 

T(z)\m 3 (z)= 16V2 16 (63) 
V 8V2 36 / 

and the three algebraic eigenfunctions have been plotted in Figure 3. It is im- 
mediately apparent that the three algebraic functions have 0, 1 and 3 nodes, 
and therefore one eigenfunction of the potential (the one with two nodes) 
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Figure 3: The three algebraic eigenfunctions ipo, 4>i and ip2 of Us(x). 



is missing from the algebraic sector, i.e. it cannot have the factorized form 
(JSOj). The situation is repeated for higher values of n and our calculations 
suggest that in the general case, the eigenfunction with n — 1 nodes does not 
belong to the algebraic sector. It remains an open question to find out what 
form this "missing" eigenfunction might have. This is the first example of a 
quasi-exactly solvable potential in which the algebraic eigenfunctions are not 
the lowest lying in the energy spectrum. 

5 Summary and conclusions 

We have shown how a more direct approach can be used to construct the 
most general differential operator of a given order r that leaves invariant a 
polynomial subspace generated by monomials. This direct approach can be 
used to address other problems that would be very difficult to treat in the Lie- 
algebraic approach, such as finding the most general differential operator that 
will map a polynomial module V n into a proper subset V n -k C V n . Problems 
of this kind will be treated in our forthcoming publication [9], where this 
approach is also extended to deal with multivariate polynomials and non- 
linear differential operators. The next step in difficulty involves extending 
the analysis to non-monomial modules (i.e. those which do not have a basis 
of monomials), a problem that remains still open. The applications of the 
direct approach to quasi-exact solvability in Quantum Mechanics point in 
two directions. First, they provide a very natural explanation of why some 
hamiltonians posses multiple algebraic sectors. As an example, we discuss 
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a generalization of the Lame potential (J38|) . which could have interesting 
applications in solid state physics. Second, we emphasize that some quasi- 
exactly solvable potentials exist which are not Lie-algebraic, as illustrated 
by the modified QES Morse potential (JoT?|) . In this example there is always 
one bound state missing from the algebraic sector, which constitutes a new 
phenomenon in the theory of quasi-exactly solvable problems. 
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